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QUASI-EXACTLY-SOLVABLE DIFFERENTIAL EQUATIONS 



ALEXANDER TURBINER 



Abstract. A general classification of linear differential and finite-difference 
operators possessing a finite-dimensional invariant subspace with a polyno- 
mial basis is given. The main result is that any operator with the above 
property must have a representation as a polynomial element of the universal 
enveloping algebra of the algebra of differential (difference) operators in finite- 
dimensional representation. In one-dimensional case a classification is given 
by algebras si2(R) (for differential operators in R) and sl'2(R)q (for finite- 
difference operators in R), osp(2, 2) (operators in one real and one Grassmann 
variable, or cquivalently, 2x2 matrix operators in R) and gl2 (R)jf ( for the 
operators containing the differential operators and the parity operator). A 
classification of linear operators possessing infinitely many finite-dimensional 
invariant subspaces with a basis in polynomials is presented. 
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ALEXANDER TURBINER 



This Chapter is devoted to a description of a new connection between Lie algebras 
and linear differential equations. 

The main idea is surprisingly easy. Let us consider a certain set of differential 
operators of the first order 



J a {x) = a a ^(x)d ll + b a (x) , 9^=-^ 



(0) 



a = 1, 2, . . . , k , x G R n , n — 1,2, ... ,n and a a '^(x), b a {x) are certain functions 
on R n . Then assume that the operators form a basis of some Lie algebra g of the 
dimension k — dimg. Now take a polynomial h in generators J a (x) and ask a 
question: 

Does the differential operator h(J a (x)) have some specific properties, 
which distinguish this operator from a general linear differential opera- 
tor? 

Generically, an answer is negative. However, if the algebra g is taken in a finite- 
dimensional representation, the answer becomes positive : 

The differential operator h(J a {x)) does possess a finite-dimensional in- 
variant subspace coinciding a representation space of the finite-dimensional 
representation of the algebra g of differential operators of the first or- 
der. If a basis of this finite-dimensional representation space can be 
constructed explicitly, the operator h can be presented in the explicit 
block-diagonal form. 

Such differential operators having a finite-dimensional invariant subspace with an 
explicit basis in functions can be named quasi-exactly-solvable. 

Up to our knowledge and understanding th e first expl icit examples of quasi- 
exactly-solvabl e pro blems were found by Razavy [ 1980 , 1981 and by Singh- Rampal- 
Biswas-Datta [ 198C ] . In an explicit form a general i dea of quasi-ex actly-solvability 
had been formulated at the first time in Turbiner [ 1988a , 1988b] and it had led 
to a complete catalogue of one-dimensional, quasi-exactly-sol vable S chroedinger 
operators in connection with spaces of polynomials Turbiner |l988c |. The term 
"quasi-exactly-solvability" has been suggested in Turbiner- Ushveridze [1987]. A 
connection of quasi-exactly-solvability and finite-dimensional representations of sfo 
was mentioned at the first time by Zaslavskii-Ulyanov [ 1984 1 . Later, the idea of 
quasi-exactly-solvability was generalized to multidimensional differential operators, 
matrix differential operators (Shifman and Turbiner [1989]), finite-difference opera- 
tors (Ogievetsky and Turbiner [ 1991 1 ) and, recently, to "mixed" operators contain- 
ing differential operators and permutation operators (Turbiner [ 1994a |). 

In Morozov et al [ 1990 1 it was described a connection of quasi-exactly-solvability 
with conformal quantum field theories (see also Halpern-Kiritsis 1989| ), while re- 
cent development can be found in Tseytlin [1994] and also in a brief review done 
by Shifman [1994|. Relationship with solid-state physics is given in (Ulyanov and 
Zaslavskii 1992| ), while recently very exciting results in this direction were found 
by Wiegmann-Zabrodin [1994a, 1994b |. A general survey of a phenomenon of the 
quasi-exactly-solvability was done in Turbiner [19941]. 



Given Chapter will be devoted mainly to a description of quasi-exactly-solvable 
operators acting on functions in one real (complex ) variable. 
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1. Generalities 

Let us take n linearly-independent functions fi(x), fa(x) . . . f n (x) and form a 
linear space 

?n(x) = (h{x)J 2 {x),...J n {x)), (1.1) 

where n is a non-negative integer and x G R. 

Definition 1.1. Two spaces Tn\x) and Tn){x) are named equivalent spaces, if 
one space can be obtained through another one via a change of the variable and/or 
the multiplication on a some function (making a gauge transformation), 

^(x)=g(x)^(y(x)) (1.2) 

Choosing a certain space (1.1) and considering all possible changes of the vari- 
able, y(x) and gauge functions, g(x), one can describe a whole class of equivalent 
spaces. It allows to introduce a certain standartization: fi(x) = 1, f2(x) = x, since 
in each class of spaces, one can find a representative satisfying the standartization. 
Thus, hereafter the only spaces of the following form 

T n {x) = (l,x,/i(x),/ 2 (x),..., /„(«)), (1.3) 

are taking into account. 

It is easy to see that once an operator h(y, d y ) acts on a space !Fn\y), one can 
construct an operator acting on an equivalent space Tn \y{x)) 

Kv, d y) = 9(x)h{x 1 d x )g~ 1 (x)\ x=x{y) , d x = . (1.4) 

Later it will be considered one of the most important particular case of the space 
(1.3): the space of polynomials of finite order 

V n+1 (x) = (l,x,x 2 ,...,x n ), (1.5) 

where n is a non-negative integer and x € R. It is worth noting the important 
property of an invariance of the space (1.5): 

x n V n+1 {-) = V n+1 {x), (1.6) 
x 

It stems from an evident feature of polynomials: let p n {x) € V n +i(x) is any poly- 
nomial, then x n p n (l/x) remains a polynomial with inverse order of the coefficients 
then in p n (x). 

2. Ordinary differential equations 
2.1. General consideration. 

Definition 2.1. Let us name a linear differential operator of the fcth order, Tf.(x, d x ) 
quasi- exactly- solvable, if it preserves the space V n +\- Correspondingly, the operator 
Ek{x, d x ) is named exactly- solvable, if it preserves the infinite flag V\ C Vi C V3 C 
■ • • C V n C . . . of spaces of all polynomials: Ek(x, d x ) : Vj 1— > Vj , j = 0, 1, 



Lemma 2.1. (Turbiner [1994b|J (i) Suppose n > (k — 1). Any quasi- exactly- 



solvable operator Tk can be represented by a k-th degree polynomial of the operators 

Jin ^ ^'X " 5 
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J° =xd x -\, (2.1) 

J'n ~ dx , 

(the operators (2.1) obey the sl2(R) commutation relations: [J^ , J°] = ±J^, 
\J + ,J~] = — 2J° FL). If n < (k — I), the part of the quasi- exactly- solvable oper- 
ator Tk containing derivatives up to order n can be represented by an nth degree 
polynomial in the generators (2.1). 

(ii) Conversely, any polynomial in (2.1) is quasi- exactly solvable. 

(Hi) Among quasi- exactly- solvable operators there exist exactly- solvable operators 
Ek C Tfc. 

Definition 2.2. Let us name a universal enveloping algebra U g of a Lie algebra 
g the algebra of all ordered polynomials in generators J ±! °. The notion ordering 
means that in any monomial in generators J ±,() all J + are placed to the left and 
all J~ to the right. 

Comment 2. 1 . Notion the universal enveloping algebra allows to make a state- 
ment that Tfc at k < n + 1 is simply an element of the universal enveloping algebra 
U s i 2 (r) of the algebra s/ 2 (R) taken in representation (2). If k > n + 1, then TJ. is 

represented as an element of J7 s ; 2 (r) plus B j x * +i , where B is any linear differential 
operator of order not higher than (k — n — 1). In otherwords, the algebra of dif- 
ferential operators acting on the space (1.4) coincides to the universal enveloping 
algebra C/ s ; 2 (r) of the algebra shCR) taken in representation (2.1) plus operators 
annihilating (1.4). 

Comment 2.2 . The algebra (2.1) has the following invariance property 

x- n Jt°(x,d x ) x n | x= i J±<°(z,d z ) 

as a consequence of the invariance (1.6) of the space V n +i- In particular, if z = 1 Jx, 
then 

J+(x,d x ) => -J-(z,d z ) 
J%(x,d x )^-J%(z,d z ) 

J-(x,d x )^-J+(z,d z ) (2.2) 

Let us introduce the grading of generators (2.1) as follows. It is easy to check 
that any s^2(R)-generator maps a monomial into monomial, J™s p i— > x p+da . 

Definition 2.3. The number d a is named a grading of the generator J" : deg(J") = 
d a . 

Following this definition 

deg(J+) = +1 , deg(J°) = , deg(J-) = -1, (2.3) 

and 

deg[{J+) n +(J° n r{Jn) n -] = n + -n_. (2.4) 

Notion of the grading allows us to classify the operators in the Lie-algebraic 
sense. 



lr T he representation (2.1) is one of the 'projectivized' representations (see Turbiner [1988a, 
1988b]). This realization of sZ2(R) has been derived at the first time by Sophus Lie. 
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Lemma 2.2. A quasi-exactly-solvable operator Tk C t^ s j 2 (R,) has no terms of posi- 
tive grading, if and only if it is an exactly-solvable operator. 

Comment 2.3 . Any exactly-solvable operator having term of negative grading pos- 
sesses terms of positive grading after transformation (2.2). A quasi-exactly-solvable 
operator always possesses terms of positive grading as in s-space representation, as 
in z-space representation. 



Theorem 2.1. (Turbiner [19941:],) Letn be a non-negative integer. Take the eigen- 
value problem for a linear differential operator of the kth order in one variable 

T k (x,d x )cp(x) = eip(x) , (2.5) 

where Tk is symmetric. The problem (2.5) has (n + 1) linearly independent eigen- 
functions in the form of a polynomial in variable x of order not higher than n, if 
and only if Tk is quasi-exactly-solvable. The problem (2.5) has an infinite sequence 
of polynomial eigenf unctions, if and only if the operator is exactly- solvable. 

Comment 2.4 ■ The " if " part of the first statement is obvious. The " only if " 
part is a direct corollary of Lemma 2.1 . 

This theorem gives a general classification of differential equations 

fe 

^2aj(x)dicp{x) = scp(x) (2.6) 

3=0 

having at least one polynomial solution in x. The coefficient functions aj{x) must 
have the form 

k+j 

i=0 

The explicit expressions (2.7) for coefficient function in (2.6) are obtained by the 
substitution (2.1) into a general, kth degree polynomial element of the universal 
enveloping algebra U s [ 2 ^y Thus the coefficients Oj^ can be expressed through the 
coefficients of the fcth degree polynomial element of the universal enveloping algebra 
^si 2 (R)- The number of free parameters of the polynomial solutions is defined by 
the number of parameters characterizing a general fcth degree polynomial element 
of the universal enveloping algebra U s i 2 rjiy A rather straightforward calculation 
leads to the following formula 

par(T k ) = (fc + 1) 2 (2.8) 

where we denote the number of free parameters of operator Tk by the symbol 
par(Tk). For the case of an infinite sequence of polynomial solutions expression 
(2.7) simplifies to 

3 

a ji x ) = X]^,^ ( 2 - 9 ) 

i=0 

in agreement with the results by Krall [ |1938| (see also Littlcjohn | 1988[ ). In this 
case the number of free parameters is equal to 

, x (fc + l)(fc + 2) 
par(Ek) = ± >j '- (2.10) 

On can show that the operators Tk with the coefficients (2.9) preserve a finite flag 
Vq C V\ C V-2 C • • • C Vk of spaces of polynomials. One can easily verify that the 
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preservation of such a finite flag of spaces of polynomial leads to the preservation 
of an infinite flag of such spaces. 

A class of spaces equivalent to the space of polynomials (1.5) is presented by 

(a(z),a(z)(3(z),...,a(z)f3(z) n ) , (2.11) 

where a(z),(3(z) are any functions. Linear differential operators acting in the 
(2.11) are easily obtained from the quasi-exactly-solvable operators (2.5)-(2.7) (see 
Lemma 2.1) making the change of variable x — (3(z) and the "gauge" transforma- 
tion T = a(z)Ta(z) _1 and have the form 

k k+j 

T k = o( 2 )J]Ea jlj ^) i )(4)| I=/SW a( Z )- 1 (2.12) 

j=0 i=0 

where the coefficients aji are the same as in (2.7). 

So the expression (2.12) gives a general form of the linear differential operator of 
the fcth order acting a space equivalent to (1.5). Since any one- or two-dimensional 
invariant sub-space can be presented in the form (2.11) and can be reduced to (1.5), 
the general statement takes place: 

Theorem 2.2. There are no linear operators possessing one- or two-dimensional 
invariant sub-space with an explicit basis other than given by Lemma 2.1. 

Therefore an eigenvalue problem (2.5), for which one eigenfunction can be found 
in an explicit form, is related to the operator 

T^ = B(x,d x )d x + q (2.13) 

and its all modifications (1.4), occurred after a change of the variable and a "gauge" 
transformation. A general differential operator possessing two eigenfunctions in an 
explicit form is presented as 

T< 2 > = B{x,d x )d 2 x + q 2 (x)d x + qi (x), (2.14) 

plus its all modifications owing to a change of the variable and a "gauge" trans- 
formation. Here B(x,d x ) is any linear differential operator, qo € R and qi, 2 {x) 
are the first- and second-order polynomials, respectively, with coefficients such that 
<l2{x)d x + qi(x) can be expressed as a linear combination of the generators J^'° 
(see (2.1)). 

2.2. Second-order differential equations. The second-order differential equa- 
tions play exceptionally important role in applications. Therefore, let us consider in 
details the second-order differential equation (2.5) possessing polynomial solutions. 
From Theorem 2.1 it follows that the corresponding differential operator must be 
quasi-exactly-solvable and can be represented as 

T 2 = C++J+J+ + c +a J+J° + c + -J+J n + c - J n J n + c — J n J n + 

c+J+ + coJ° + C-J~ +c, (2.15.1) 

where c Q /3,c Q ,c e R. The number of free parameters is par(T 2 ) — 9. Under the 
condition c ++ = c +0 = c + = 0, the operator T 2 becomes exactly- solvable (see 
Lemma 2.2) 

E 2 = C+-J+J- + Co-J®J~ + C—J-J- + c J° + c-J- + c, 

(2.15.2) 
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and the number of free parameters is reduced to par(E 2 ) = 6. 

Lemma 2.3. If the operator (2.15.1) is such that 

n 

C++ = and c+ = (— — m)c+o , at some m = 0, 1, 2, . . . 

2 (2.16) 

then the operator T 2 preserves both V n +i and V m +i- In this case the number of 
free parameters is par(T 2 ) = 7. 

In fact, Lemma 2.3 claims that T 2 (J", c a p, c a ) can be rewritten as T 2 (J^, c' Q(3 , c' a ). 
As a consequence of Lemma 2.3 and Theorem 2.1, in general, among polynomial 
solutions of (2.6) there are polynomials of order n and order m. 

Remark. From the Lie-algebraic point of view Lemma 2.3 means the existence 
of representations of second-degree polynomials in the generators (2.1) possessing 
two invariant sub-spaces. In general, if n in (2.1) is a non-negative integer, then 
among representations of fcth degree polynomials in the generators (2.1), lying in 
the universal enveloping algebra, there exist representations possessing 1,2, ...,k 
invariant sub-spaces. Even starting from an infinite-dimensional representation of 
the original algebra (n in (2.1) is not a non-negative integer), one can construct the 
elements of the universal enveloping algebra having finite-dimensional representa- 
tion (e.g., the parameter n in (2.16) is non-integer, however T 2 has the invariant 
sub-space of dimension (m + 1))- Also this property implies the existence of repre- 
sentations of the polynomial elements of the universal enveloping algebra ?/ s ; 2 (r), 
which can be obtained starting from different representations of the original algebra 
(2.1). 

Substituting (2.1) into (2.15.1) and then into (2.6), we obtain 

- P A (x)c£<p(x) + P 3 (x)d x <p(x) + Pi{x)ip{x) = e<p{x), (2.17) 

where the Pj(x) are polynomials of jth order with coefficients related to c a p,c a 
and n. In general, problem (2.17) has (n + 1) polynomial solutions of the form of 
polynomials in x of nth degree. If n = 1, as a consequence of Theorem 2.2, a more 
general spectral problem than (2.17) arises (cf. (2.14)) 

- F 3 (x)d 2 x tp(x) + q 2 (x)d x ip(x) + qi(x)ip(x) = eip(x), (2.18) 

possessing only two polynomial solutions of the form (ax + b). Here F 3 is an 
arbitrary complex function of x and qj(x),j = 1,2 are polynomials of order j the 
same as in (2.14). For the case n — (one polynomial solution, p = const) the 
spectral problem (2.5) becomes (cf. (2.13)) 

- F 2 (x)d 2 x tf(x) + F 1 (x)d x <p(x) + qo<p(x) = e<p(x), (2.19) 

where F 2t i(x) are arbitrary complex functions and go G R- 

Substituting (2.1) into (2.15.2) and then into (2.6), we obtain 

- Q 2 (x)dl<p(x) + Qi(x)d x tp(x) + Qo(x)ip(x) = e<p(x), (2.20) 

where the Qj(x) are polynomials of jth order with arbitrary coefficients. One can 
easily show that the differential operator in the r.h.s. can be always presented in 
the form (2.15.2). The coefficients of Qj(x) are unambiguously related to c a p,c a 
for any value of the parameter n. Thus n is a fictitious parameter and, for instance, 
it can be put equal to zero. 
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2.3. Quasi-exactly-solvable Schroedinger equations (examples). From the 
point of applications the Schroedinger equation 

(-4 + V{z))9(z) = e*(z), (2.21) 

is one of the important among the second-order differential equations (see e.g. 



Landau and Lifschitz [1974Q. Here e is the spectral parameter and ^(z) must 
be square-integrable function on some space. It can be the whole real line, semi- 
line or a finite interval. Therefore it is quite natural to search the quasi-exactly- 
solvable and exactly-solvable operators of the Schroedinger type acting on a finite- 
dimensional linear space of square-integrable functions:/i(z), i = 1,2, . . .n. 

One possible way to get the quasi-exactly-solvable and exactly-solvable Schroedinger 
operators is to transform the quasi-exactly-solvable Ti and exactly-solvable E 2 oper- 
ators acting on finite-dimensional spaces of polynomials (1.5) into the Schroedinger 
type operators. It always can be done making a change of a variable and a gauge 
transformation (sec (1.2)) as a consequence of one-dimensional nature of equations 
studied. In practice, the realization of this transformation is nothing but a conver- 
sion of (2.17)-(2.20) into (2.21). The only open question remains: does new basis 
belong to square-integrable one or not? This question will be discussed in the end 
of this Section. In the following consideration we restrict ourself by the case of real 
functions and real variables. 

Introducing a new function 

= v{x(z))e- A{z \ (2.22) 

and new variable x — x(z) , where A(z) is a certain real function, one can reduce 
(2.17)-(2.20) to the Sturm-Liouville-type problem (2.21) with the potential equal 
to 

V{z) = (A') 2 -A" + P 2 {x{z)) , (2.23) 

if 

,P%, , , , f dx 



A = ( — )dx — logz' , z = ± . 



for the case of (2.17), or 



A = I (^-)dx — logz' , z — ± / —== 



for the case of (2.18), or 



A = I (— -r)dx — logz 1 , z = ± ' 



iV " ' J VF2 

for the case of (2.19) with replacement of P2(x(z)) for two latter cases by Q±(x(z)) 
or Q , respectively. If the functions (2.22), obtained after transformation, belong 
to the £ 2 (2? )-space[] , we arrive at the recently discovered quasi-exactly-solvable 



Schroedinger equations (Turbincr [1988a, 1988b, 1988c|), where a finite number of 
eigenstates is found algebraically. 

In order to proceed a description of concrete examples of quasi-exactly-solvable 
Schroedinger equations, first of all let us generalize the eigenvalue problem (2.21) 
inserting a weight function g(z) = g(x{z)) in the r.h.s. 

(-4 + V(z))9(z) = eg(x(z))^(z), (2.24) 

2 Depending on the change of variable x = x(z), the space X> can be whole real line, semi-line 
and a finite interval. 
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where for a sake of future convenience the weight function is presented in the form 
of composition g(x(z)). This equation can obtained from (2.17)-(2.19) by taking 
the same gauge factor (2.22) as before but with another change of the variable 

dx 



s /g(x)P i (x) 

what leads to a slightly modified potential then in (2.23) 

V(z) = (A') 2 - A" + P 2 (x(z))q(x(z)) 

Below we will follow the catalogue given at Turbiner ]1988c( |. A presentation of 
results is the following: firstly, we display the quadratic element T 2 of the universal 
enveloping algebra sl 2 in the representation (2.1) and its equivalent form of differen- 
tial operator T 2 (x, d 2 ), secondly, the corresponding potential V(z)a,nd afterwards, 
the explicit expression for the change of the variable x — x(z), the weight function 
g(z)and finally the functional form of the eigenfunctions ^(z) of the "algebraized" 
part of the spectra 0. 

We begin a consideration with the quasi-exactly-solvable equations, associated 
to the exactly-solvable Morse oscillator; it implies that at the limit, when number 
of "algebraized" eigenstates (n + 1) goes to infinity, the Morse oscillator occurs. 

Comment 2.5 The Morse oscillator is one of well-known exactly-solvable quantum- 



mechanical problems (see e.g. Landau and Lifschitz [ 1974 1 ) . It is described by the 
Schroedinger operator with the potential 

V{z) = A 2 e~ 2az - 2Ae- az , A, a > . 

This potential is used to model the interaction of the atoms in diatomic molecules. 



T 2 = -a 2 J+J~ + 2aaJ+ - a[a{n + 1) + 2b}J° - 2acJ~ (2.25) 
- — [a(n + l) + 2b] 

or as the differential operator, 

T 2 (x, d x ) — —a 2 x 2 d 2 + a[2ax 2 — (a + 2b)x — 2c]d x — 2aanx 

leads to 

V(z) = a 2 e- 2az - a[2b + a(2n + l)] e - az + c{2b - a)e az + c 2 e 2az 



+b 2 - 2ac 

where 

x = er az , g = 1 , 
with the eigenfunctions of "algebraized" part of the spectra 

9( z ) = Pn ( e ~ az )exp(--e- az - bz - -e az ) 

a a 

at a > 0, a, c > and V6. 



(2.26) 



3 The functions p n (x) occurring in the forthcoming expressions for ^(z) are polynomials of the 
nth order. They are nothing but the polynomial eigenfunctions of the operator T2(x,d x ) . 
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II. 



T 2 = -aJ° n J- + 2aJ+ - 2cJ° - [a(- + 1) + 26] J" - cn (2.27) 

or as the differential operator, 

T-zix, d x ) = —axd 2 + (—2ax 2 + 2cx + a + 2b)d x + 2anx 

leads to 

V(z) = a 2 e- 4az - 2ace- 3az + [c 2 - 2a{b + aan + a)] e - 2az + c{2b + a)e~ az + b 2 , 



(2.28) 



where 



x = e- az , q= -e~ az 
a 



= p n (e- az ) exp ( e - 2az + -e- az - bz) 

2a a 

at a > 0, a > 0, 6 > and Vc. 
///. 

T 2 = -aJ+J° + (26 - 3an/2) J+ - 2aJ° - 2cJ- - an (2.29) 
or as the differential operator, 

T 2 {x, d x ) — —ax 3 d 2 + [(2b — a)x 2 — 2ax — 2c]d x + (an — 2b)nx 

leads to 

V(z) = c 2 e iaz + 2ace 3az + [a 2 - 2c(b + a)]e 2az - a(2b + a)e az + b 2 + an(an - 2b) , 

(2.30) 

where 

x = e~ az , q = -e az , 
a 

= p n (e- az )exp(-^e 2az - -e az + bz) 
2a a 

at a > 0, b > 0, c > and Va. 

Next two quasi-exactly-solvable potentials are associated to the one-soliton or 
Poschle- Teller potential. 

Comment 2. 6 The Poschle- Teller or one-soliton potential describes another well- 
known exactly-solvable quantum-mechanical problem (see e.g. Landau and Lifs- 



chitz 1974 ). It is given by the Schroedinger operator with the potential 

1 
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V(z) = -A*— — . 

cosn az 

This potential has a unique property of an absence of reflection and is also one of 
the simplest solutions of so called Korteweg-de Vries equation playing important 
role in the inverse problem method (for detailed discussion see e.g. the book by 
V.E. Zakharov et al |1980Q ). 
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IV. 

T 2 = -Aa 2 J+,f n + 4a 2 J+J- - 2a[a(3n + 2k + 1) + 2c] J+ 
+ 2a[a(n + 2) + 2c - 2a] J° + 4aa,J~ + an[a(n + 2) + 2c - 2a] 

(2.31) 

or as the differential operator, 

T 2 (x, d x ) = -4a 2 (x 3 - x 2 )d\ - 2a[(3a + 2ak + 2c)x 2 - 2(a - a + c)x - 2a]d x 
+ (2n + k)a[{2n + k + I) a + 2c]x 

leads to 

V(z) = a 2 cosh 4 az — a(a + 2a — 2c) cosh 2 az (2.32) 
— [c(c + a) + a(2n + k)(a(2n + k) + a + 2c)] cosh -2 az + c 2 + aa — 2ac 

where 

x = cosh az -2 ,£>=1, 

^(z) — (tanhaz) fe p n (tanh 2 az) (cosh az) _c / Q exp (— — cosh2az) 

4a 

at a > 0, a > 0, Vc and k = 0, 1. 
V. 

T 2 = -Aa 2 J+J- + 4a 2 J°J- + 4abJ+ - 2a[a{2n + 2k + 3) + 2a + 46] J° 
+ 2a[a(n + 2) + 2a + 26) J" - a[an(2n + 2k + 3) + 2an - 2bk] 

(2.33) 

or as the differential operator, 

T 2 (x, dx) = -Aa 2 x{x - l)d 2 x + 2a[2bx 2 - (2a + 46 + 2ka + 3a)x + 2(a + a + b)]d x 

-Aabnx + 2ab(2n + k) 

leads to 

V(z) = -b 2 cosh -6 az + b[2a + 36 + a(4n + 2k + 3)] cosh -4 az 

(2.34) 

-[(a + 36) (a + 6 + a) + 2(2n + k)ab)] cosh -2 az + (a + b) 2 

where 

x = cosh az , q = cosh az , 

, o -(a + b) 6 9 

^(z) = (tanhaz) p„(tanh az)(coshaz) = exp (— tanh 2az) 
at a > 0, (a + 6) > and k = 0, 1. 

Next two quasi-exactly-solvable potentials are associated to the harmonic oscil- 
lator potential. 

VI. 



12 



ALEXANDER TURBINER 



It is the first example of the quasi-exactly-solvable Schroedinger operator. Let us 



take the following non-linear combination in the generators (2.1) (Turbiner [1988c | ) 
T 2 = -4J° J,7 + 4aJ+ + AbJ° - 2(n + 1 + 2k) J- + 2bn (2.35) 

or as the differential operator, 

T 2 (x, d x ) = -Axd 2 x + 2(2ax 2 + 2bx - 1 - 2k)d x - Aanx , 

where x G R and a > 0, V& or a > 0, b > . Putting x = z 2 and choosing the gauge 
phase A = ax 2 /4 + bx/2 — /c/21na;, we arrive at the the spectral problem (2.19) 



with the potential (Turbiner and Ushveridze 1 1987 ] ) 



V{z) = a 2 z b + 2abz 4 + [V - (4n + 3 + 2k)a]z 2 - 6(1 + 2k), (2.36) 

for which at k = (k = 1) the first (n + 1) eigenfunctions, even (odd) in x, can be 
found algebraically. Of course, the number of those "algebraized" eigenfunctions 
is nothing but the dimension of the irreducible representation (1.4) of the algebra 
(2.1). {n + 1) 'algebraized' eigenfunctions of (2.19) have the form 

= zV(z 2 Kt"t ] (2.37) 

where p n (y) is a polynomial of the nth degree. It is worth noting that if the pa- 
rameter a goes to 0, the potential (2.35) becomes the harmonic oscillator potential 
and polynomials z k p n (z 2 ) reconcile to the Hermite polynomials H2 n +k{z) (see dis- 
cussion below). 

VII. 



T 2 = -47° J' + AaJ+ + 46J° - 2(n + d + 21 - 2c) J- + 2bn 

(2.38) 

or as the differential operator, 

T 2 (x, d x ) = -Axd 2 x + 2(2ax 2 + 2bx - d - 21 + 2c)d x - Aanx 

leads to 

V(z) = a 2 z 6 + 2abz 4 + [b 2 - {An + 21 + d + 2 - 2c)a]z 2 + c(c - 21 - d + 2)z~ 2 

(2.39) 

-b(d + 2l-2c), 



where 



Z 2 , Q=l 



*(z) = p n (z 2 )z l C C *T- 

at a > 0, V6 or a > 0, b > and (d + 1 — c) > 1 . This case corresponds to the radial 
part of (i-dimensional Schroedinger equation, z G [0,oo). At d — 1 the radial part 
coincides with the ordinary Schroedinger equation (2.19) at z G R. The potential 
(2.39) becomes a generalized version of the potential (2.36) with an additional 
singular term proportional to z~ 2 . 



Next two quasi-exactly-solvable potentials are associated to the Coulomb prob- 
lem. 



LIE-ALGEBRAS AND QUASI-EXACTLY 



13 



VIII. 



T 2 = -J°Jn + ZaJn + 2bJ n ~ O/ 2 + d + 21 - 2c - 1) J" - an 

(2.40) 

or as the differential operator, 

T 2 (x, d x ) = -xd% + {2ax 2 + 2bx + 2c - d - 21 + l)d x - 2anx 

leads to 

V(z) = a 2 z 2 + 2abz - 6(2/ + d - 1 - 2c)z~ 1 + c(c - 21 - d + 2)z" 2 

(2.41) 

+b 2 + a(2c - d - 21 - 2n) , 

where 

= p„(z)z^ c e-^ 2 - b2 , 
at a > 0, 6 > and (ef + / — c) > 1. 

IX 



T 2 = -J+Jn + 2aJ+ -(n + d-1 + 21- 2c)J° + 26J- - n(d + 21 - 1 - 2c) 

(2.42) 

or as the differential operator, 

T 2 (x, d x ) = -x 2 d\ - [2ax 2 + (2c-d-2l + l)x - 2b]d x - Aanx 

leads to 

V[z) = 6 2 z~ 4 - b(2c -2l~d + 3)z~ 3 + [c(c - 21 - d + 2) - 2afe]z" 2 

(2.43) 

-a(2n + 2l + d-l- 2c)z~ 1 + a 2 , 

where 

= p n (z)z l - c e- az ~ hz ~\ 
at a > 0, b > or a > 0, b > and c,l,d£ R. 

Now let us show an example of the non-singular periodic quasi-exactly-solvable 
potential connected with Mathieu potential. 

Comment 2.7 The Mathieu potential 

V(z) = Acosaz 

is one of the most important potentials in many branches of physics, engineering etc. 
Detailed description of the properties of the corresponding Scroedinger equation 



can be found, for instance, in Bateman-Erdelyi [1953], vol.3 and also Kamke [ 1959 1 , 
Equation 2.22. 
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X. 

Firstly, take the following operator 

T2 = ^ J+_ M J"_ M - a 2 J-_^ J-_^+ (2.44) 

2aaJ+_ M + a 2 (n + /i + 1) J°_„ - 2 M J„^ + a 2 ^ ^ ^ 

or as the differential operator, 

^2(25, d x ) = a 2 (x 2 — l)d 2 + a[2ax 2 + a(l + 2/i)x — 2a]d x — 2aa(n — fx)x 
The transformation (2.22)-(2.23) leads to the periodic potential 

V(z) = a 2 sin 2 az - aa(2n + 1) cos az + fia 2 , ( 2 -45) 

where 

x = cos az , q = 1 , 
*0) = (sinaz)>„_ p (cosaz) e ( = COSQz) , 
at a > 0, a > 0. Here fi — 0,1 . For the fixed n, (2n + 1) eigenstates having a 
meaning of the edges of bands can be found algebraically. 
Secondly, take the following operator 

T 2 = otJ+^J-^ - a 2 J-_ x J-_^ (2.46) 

2 (n 2 -l) 

a 



2 



2aaJ^_ x + a 2 (n + 1) J„_ t — a[2a + a{v x — vz)]J^-i 
or as the differential operator, 

T 2 (x, dx) — a 2 (x 2 — l)d 2 + a[2aa; 2 + 2ax — 2a — <x{y\ — i/ 2 )]d x — 2aa(n — l)x 
The transformation (2.22)-(2.23) leads to the periodic potential 

a 2 

V(z) = a 2 sin 2 az — aa(2n) cos az + aa(v\ — 1/2)] — , (2.47) 

where 

x = cos az , q = 1 , 

$(z) = (cosaz) ,yi (sinaz) ,y2 p„_i(cosaz) e ^ cosaz \ 

at a > 0, a > 0. Here 1^,2 = 0, 1, but 1/1 + 1/3 = 1. For the fixed n, (2n + 1) 
eigenstates having a meaning of the edges of bands can be found algebraically. 

2.4. Quasi-exactly-solvable Schroedinger equations (Lame equation). In 

this Section we consider one of the most important second-order ordinary differen- 
tial equation - m-zone Lame equation 

-d 2 x ^ + m(m+ l)V(x)^ = (2.48) 

where V{x) is the Weierstrass function in a standard notation (see, e.g. Bateman 
and A. Erdelyi [ 1953 ) , which depends on two free par amete rs, and m = 1, 2, . . . . 
In a description we mainly follow the paper (Turbiner | 1989 |). 

The Weierstrass function is a double-periodic meromorphic function for which 
the equation V' 2 = (V — ei)(V — e 2 )(V — ^z) holds, where ^ e,; = 0. I ntroducing 
the new variable £ = Vix) + | X) a * m (( 2-48 )) (see, e.g. Kamke [1959Q, the new 
equation emerges 

,.1,1 1 1 . . m(m + l)£ + e 

V + 9 7 + 7 + 7 h - 777 T77 \77 = 

2 S-ffli ?-a2 ?-a3 4(^-a 1 )(^-a 2 )(^-a 3 ) 
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where rj(£) = ip( x )- Here the new parameters a, satisfy the system of linear equa- 
tions 6j = aj — i ^2 a.i . Equation ( |2.49 ) is named by an algebraic form for the Lame 
equation. There exists a spectral parameter A for which equation (2.49) has four 
types of solutions: 



4 2) = {Z-aifftpk® , i = 1,2,3 
4 3) =K-oi 1 1/2 (e-ai») 1/2 Pfc-i(O. h¥=l2;i^h,2;i,h,2 



(2.50.1) 

(2.50.2) 

1,2,3 
(2.50.3) 

(2.50.4) 

where Pr(£) are polynomial in £ of degree r. If the value of parameter n is fixed, 
there are (2m + 1) linear independent solutions of the following form: if m = 2k 
is even, then the j/ 1 ^) and r]( 3 \£) solutions arise, if m = 2fc + 1 is odd we have 
solutions of the j/ 2 ^(£) and f/ 4 ^(£) t ypes . Those eigenvalues have a meaning of the 
edges of the zones in the potential fl2.< 



(4) = K - «i) 1/2 (£ - a 2 ) 1/2 (e - oaJpfc-xCO 



Theorem 2.3. (Turbiner 1 1989 1) Tfee sp ec tral pro ble m (2.4& ) at m — 1,2,... 
with polynomial solutions ( 2.50.1 ), ( 2.50.^ ), \2.50.S ), ( 2.50.4 ) * s equivalent to the 
spectral problem (2.5) for the operator Ti (2.15.1) belonging the universal enveloping 
sli-algebra in the representation (2.1) with the coefficients 

c+o = 4 , c + _ = -4^a 4 , c - = 4}^^ , c — = ai a 2 a 3 

(2.51) 

before the terms quadratic in generators and the following coefficients before linear 
in generators J^'° : 

(1) For (£) -type solutions at m — 2k, r = k 

c+ = -6k - 2 , c = 4(fc + 1) Oi > c - = -2 (A: + 1) X! a * 

(2) For 77) (£)-type solutions at m = 2fc + l,r = k 

c + = —6k — 6 , cq — A(k + 2) a% — a, , 

c_ = -2(fc + 1) ^diflj - 4a ;i a i3 , i^h,2,h^h 



(2.52.1) 



(2.52.2) 



(3) For 77) (£)-type solutions at m = 2k, r = k — 1 
c+ = 



: -6k - 4 , c = 4(fc + 1) V" a,i + 4a 4 , 
2(k + 2) ^ a i a i + 4a h ai 2 , i ^ i 1)2 , Zi 7^ Z 2 



(4) For n^^) -type solutions at m = 2k + 1, r = k — 1 



-6k- 



, c =4(fc + 2)^\ 



-2(fc + 2)£ 



(2.52.3) 



(2.52.4) 
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So, each type of solution fl2.50.lf ), J2.50.2D , d2.50.3fj , ( |2.50.4| ) corresponds to 
the particular spectral problem ( 2.4S ) with a special set of parameters ( |2.5l| ) plus 
(2.52.1), (2.52.2), (2.52.3), (2.52.4), respectively. It can be easily shown that the 
calculation of eigenvalues e of (2.48) corresponds to the solution of a characteristic 
equation for the four-diagonal matrix: 

C U -x = (l-l-2j)[4(j+l-l)+c+] , 

C u = [l(2j + l-l)c+- + (l-j)co] , 

C u+1 = (l+l)(j-l)co- + (l + l)c^ , 



C M+2 = -(/ + l)(i + 2)c 



(2.53) 



where the size of this matrix is (k + 1) x (fc + 1) and 2j = k for (2.50.1), (2.50.2), and 
k x k and 2j = k — 1 for (2.50.3), ( |2.50.4 ), respectively. In connection to Theorem 
2.3 one prove the following theorem (Turbiner [ 1989 1 ) . 



Theorem 2.4. Let fix the parameters e's fa's) in 12.48 ) (or \2.4$ )) except one, 
e.g. ei(ai). The first (2m + 1) eigenvalues of ft2-4S\ ) (or ( 2.4$ )) form (2m + 1)- 



sheeted Riemann surface in parameter ei(ai). This surface contains four discon- 
nected pieces: one of them corresponds to rf 1 ^ (n^) solutions and the others corre- 
spond to rj^\ri^). At m = 2k the Riemann subsurface for rf^ has (k + 1) sheets 
and the number of sheets in each of the others is equal to k. At m = 2k + 1 the 
number of sheets for n^ is equal to k and for rf 2 ^ each subsurface contains (k + 1) 
sheets. 

It is worth emphasizing that we cannot find a relation between the spectral 
problem for the two-zone potential 



V 



-25>(z 

fc=i 



3 

i=i 







(2.54) 



(see Dubrovin and Novikov [ 1974 1 ) ^ and the spectral problem (2.5) for T2 with the 
parameters (2.51) and (^ 52 T|) or ( |2.52.3| ) at k = 1. In this case eigenvalues e and 

also eigenfunctions (2.49) (but not (2.48) do not depend on parameters c 

Comment 2.8 . One can generalize a meaning of isospectral deformation, saying 
we want to study a variety of potentials with the first several coinciding eigenvalues. 
It can be named quasi-isospectral deformation. 



Now let us consider such a quasi-isospectral deformation of (2.48) at m = 2. 
It arises from the fact that the a ddition of t he term C++ J+ J+ to the operator T2 
with the parameters ( 2.51 ) and ( 2.52.1 ) or ( 2.52.3 ) at k = 1 does not change the 



characteristic matrix ( 2.53 ). Making the reduction (2.22)-(2.23) from the equation 



(2.5) to the Schroedinger equation (2.21), we obtain 



V(x) 



2c++£ 6 



P2(0 , 



(2.55) 



where 



4 2 + c -£ + c__ , P 2 (0 = -m(m + 1)C + tt 



co 

2(2.56) 



4 The potential (2.54) and the original Lame potential (2.48) at m = 2 are related to via the 
isospectral deformation. 
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and £ is defined via the equation 



(2.57) 



In general, the potential (2.55) contains four double poles in x and does not reduce 
to ( 2.54 ). It is worth noting that the first five eigenfunctions in the potential ( [2.55 ) 
have the form 



A£ + B 

(£-« 4 ) 1/2 (e-a,y/2 



exp 



, i¥=3, h3 = 1,2,3. 

(2.58) 



Here £ is given by (2.57). The first five eigenvalues of the potential (2.55) do not 
depend on the parameters c , c ++ . 

2.5. Exactly-solvable equations and classical orthogonal polynomials. As 

one of the most important properties of the exactly-solvable operators (2.15.2) is 
the following: the eigenvalues of a general exactly-solvable operator E 2 are given 
by the quadratic polynomial in number of eigenstate 

&m = c oo™ 2 + c o m + const (2.59) 

(for details see Turbiner fll988a| , |l988b| ). It can be easily verified by straightforward 
calculation. 

Taking different exactly-solvable operators E 2 (see (2.15.2)) for the eigenvalue 
problem (2.5) one can reproduce the equations having the Hermite, Laguerre, Le- 
gendre and Jacobi polynomials as eigenfunctions (Turbiner [1988a, 1988t]), which 
is shown below. In the definition of the about polynomials we follow the definition 
given in (Bateman and A. Erdelyi [1953]). 

1. Hermite polynomials. 

The Hermite polynomials H 2 m+k, m = 0,1,2, ... , k = 0, 1 are the polynomial 
eigenfunctions of the operator 

dl 



E 2 {x,d x ) 



2xd 

nr. 



(2.60) 



which immediately can be rewritten in terms of the generators (2.1) following the 
Lemma 2.1 



En 



2J °(z) 



(2.61) 



J (x)J (x) 

However, there exists another way to represent the operators related to the Hermite 
polynomials. Let us notice that k has a meaning of the parity of the polynomial 
H2m+k and 

H 2m +k{x) = x k h m {x 2 ) 
Then it is easy to find that the operator having h m (y) as the eigenfunctions 

E 2 (y, d y ) = Ayd 2 y - 2(2y - 1 - 2k)d y (2.62) 

and, correspondingly 

E 2 = 4J °(y)J -(y) - 4J °(y) + 2(1 + 2k) J ~ (y) (2.63) 

Of course, those two representations are equivalent, however, a quasi-exactly-solvable 
generalization can be implemented for the second representation only (see examples 
VI and VII in Section 2.3). 
2. Laguerre polynomials. 
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The associated Laguerre polynomials L 1 ^ (x) occur as the polynomial eigenfunc- 
tions of the generalized Laguerre operator 

E 2 (x,d x ) = xd 2 x + (a+1- x)d x (2.64) 



where a is any real number. Of course, the operator (2.64) can be rewritten as 

E 2 = Jo Jo - J o + (a+W (2-65) 
3. Legendre polynomials. 

The Legendre polynomials P2 m +k(x) are the the polynomial cigcnfunctions of 
the operator 

E 2 (x,d x ) = (1 - x 2 )d 2 x - 2xd x (2.66) 

or, 

E 2 = -« + Jo Jo ~ 4? (2-67) 

Analogously to the Hermite polynomials there exists another way to represent the 
operators related to the Legendre polynomials. Let us notice that k has a meaning 
of the parity of the polynomial P 2m +k and 

P2rn+k(x) = X k p m {x 2 ) 

Then it is easy to find that the operator having p m (y) as the eigenfunctions 

E 2 (y, d y ) = 42/(1 - y)d 2 y + 2[1 + 2k - (3 + 2k)y]d y (2.68) 

and, correspondent^ 

E 2 = ~4J+(y)J^(y) + 4J°(y)J^(y) - 2(3 + 2k) J °(y) + 2(1 + 2k)J (y) 

(2.69) 

4- Jacobi polynomials. 

The Jacobi polynomials appear as the polynomial cigcnfunctions of Jacobi equa- 
tion taking in cither symmetric form with the operator 

E 2 {x, d x ) = (1 - x 2 )d 2 x + [b-a-(a + b + 2)x]d x (2.70) 

corresponding to 

E 2 = - J ° J ° + Jq Jq - (1 + a + b) J ° + {b- a)J , (2.71) 



or asymmetric form (see e.g. the book by Murphy 1 1960 or Bateman-Erdelyi 



1953| ) 

E 2 (x, d x ) = x(l - x)d 2 x + [l + a-(a + b + 2)x]d x (2.72) 

corresponding to 

E 2 = - J ° J ° + J^J - (1 + a + b) J ° + (a + 1) Jq, (2.73) 

Under special choices of the general element E^Eq, Eg), one can reproduce all 
known fourth- (sixth-, eighth-)order differential equations giving rise to infinite se- 



quences of orthogonal polynomials (see e.g. Littlejohn [ 1988| and other papers in 
this volume). 



Recently, A. Gonzalez-Lopez, N. Kamran and P. Olver [1993 gave the complete 
description of second-order polynomial elements of t/ s ; 2 (R) in the representation 
(2.1) leading to the square-integrable eigenfunctions of the Sturm-Liouville problem 
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(2.21) after transformation (2.22)-(2.23). Consequently, for second-order ordinary 
differential equation (2.17) the combination of their result and Theorems 2.1, 2.2 
gives a general solution of the problem of classification of equations possessing a 
finite number of orthogonal polynomial solutions. 
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3. Finite-difference equations in one variable 
3.1. General consideration. Let us define a multiplicative finite-difference op- 



erator, or a shift operator or the so-called Jackson symbol (see e.g. Exton [ 1983 1 , 
Gasper and Rahman [ 199(| ) 

Df{x) = (l-q) X (3 ' 1} 

where q € R and f(x) is real function x S R . The Leibnitz rule for the operator 
D is 

Df{x)g{x) = (Df(x))g(x) + f(qx)Dg(x) 
Now one can easily introduce the finite-difference analogue of the algebra of the 
differential operators (2.1) based on the operator D instead of the continuous de- 



rivative (Ogievetsky and Turbiner [ 1991 1 ) 



J+ = x 2 D - {n} 



x 



J° = xD-n (3.2) 
J- = D, 



where {n} — \_H q is so called q number and h = ^r^^p ■ The operators (3.2) 
after multiplication by some factors 

, _ g- n {2n + 2} -Q 

■' p + 1 {n + 1} " 

(see Ogievetsky and Turbiner [ 1991 1 ) form a quantum algebra s?2(R) g with the 
following commutation relations 



qfj ~ j f = -f 



q 



J + r - r~j + = -(« + !)? (3.3) 



f~t q] + J° = j 



The parameter q does characterize the deformation of the commutators of the 
classical Lie algebra s?2- If <7 — ► 1, the commutation relations (3.3) reduce to the 
standard s^(R) ones. A remarkable property of generators (3.2) is that, if n is a 
non-negative integer, they form the finite- dimensional representation corresponding 
the finite- dimensional representation space V n +i the same as of the non-deformed 
sl2 (see (2.1)). For values of q others than root of unity this representation is 
irreducible. 

Comment 3.1 The algebra (3.3) is known in literature as so-called the second 



Witten quantum deformation of sl2 in the classification of C. Zachos [ 199 1 1 ) . 

Similarly as for differential operators one can introduce quasi-exactly-solvable 
Tk{x,D) and exactly-solvable finite-difference operators Ek{x,D) (see definition 
2.1). 



Lemma 3.1. (Turbiner [ 1994b ^ (i) Suppose n > (k — 1). Any quasi- exactly- 
solvable operator T^, can be represented by a kth degree polynomial of the operators 
(3.2). If n < (k — 1), the part of the quasi-exactly-solvable operator contain- 
ing derivatives up to order n can be represented by a nth degree polynomial in the 
generators (3.2). 
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(ii) Conversely, any polynomial in (3.2) is quasi- exactly solvable. 
(Hi) Among quasi-exactly-solvable operators there exist exactly-solvable operators 
EkCf k . 

Comment 3.2 If we define an analogue of the universal enveloping algebra U g for 
the quantum algebra g as an algebra of all ordered polynomials in generators, then a 
quasi-exactly-solvable operator Tj. at k < n+1 is simply an element of the 'universal 
enveloping algebra' J7 s j 2 (r,) of the algebra s?2(R) 9 taken in representation (3.2). If 
k > n + 1, then Tfc is represented as an element of £/ s z 2 (r) plus BD n+1 , where B 
is any linear difference operator of order not higher than [k — n — 1). 

Similar to sZa(R) (see definition 2.3), one can introduce the grading of generators 
(3.2) of sl2(R) q (cf. (2.3)) and, hence, the grading of monomials of the universal 
enveloping U sh{R)q (cf. (2.4)). 

Lemma 3.2. A quasi-exactly-solvable operator C ?7 s ; 2 (h.) has no terms of pos- 
itive grading, iff it is an exactly-solvable operator. 



Theorem 3.1. (Turbiner [ 119941^ Let n be a non-negative integer. Take the eigen- 



value problem for a linear difference operator of the k-th order in one variable 

f k (x,D)<p{x) = eip{x), (3.4) 

where is symmetric. The problem (3.4) has (n + 1) linearly independent eigen- 
functions in the form of a polynomial in variable x of order not higher than n, if 
and only ifT^ is quasi-exactly-solvable. The problem (3.4) has an infinite sequence 
of polynomial eigenf unctions, if and only if the operator is exactly- solvable E^. 

Comment 3.2 Saying the operator is symmetric, we imply that, considering 
the action of this operator on a space of polynomials of degree not higher than 
n, one can introduce a positively-defined scalar product, and the operator T^ is 
symmetric with respect to it. 

This theorem gives a general classification of finite-difference equations 

k 

^^aj(x)D 3 (p(x) — ecp(x) (3.5) 

3=0 

having polynomial solutions in x. The coefficient functions must have the form 

k+j 

i=0 

In particular, this form occurs after substitution (3.2) into a general fcth degree 
polynomial element of the universal enveloping algebra U s i 2 ( R ) . It guarantees the 
existence of at least a finite number of polynomial solutions. The coefficients cijj 
are related to the coefficients of the fcth degree polynomial element of the univer- 
sal enveloping algebra U s i 2 (jl) . The number of free parameters of the polynomial 
solutions is defined by the number of free parameters of a general A:-th order polyno- 
mial element of the universal enveloping algebra U s i 2 ^ q A rather straightforward 

5 For quantum sl2(R) q algebra there are no polynomial Casimir operators (see, e.g. Zachos 



1991 1). However, in the representation (3.2) the relationship between generators analogous to the 
quadratic Casimir operator 

qj+j- - J°J° + ({n + 1} - 2n) J° = h(h - {n + 1}) 
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calculation leads to the following formula 

par(f k ) = (k + l) 2 + 1 

(for the second-order finite-difference equation par(T 2 ) = 10). For the case of an 
infinite sequence of polynomial solutions the formula (3.6) simplifies to 

3 

°-]{ x ) = ^HiX* ( 3 ' 7 ) 

and the number of free parameters is given by 

, - , (k + l)(fc + 2) 
par(E k ) = ± >j >- + l 

(for k = 2, par(E 2 ) = 7). The increase in the number of free parameters com- 
pared to ordinary differential equations is due to the presence of the deformation 
parameter q. 

3.2. Second-order finite-difference exactly-solvable equations. In Turbiner 



1992 it is implemented a description in the present approach of the q-deformed 



Hermite, Laguerre, Legendre and Jacobi polynomials (for definitions of these poly- 
nomials see Exton [ 1983 1 , Gasper and Rahman | 1990[ ). In order to reproduce the 
known q-deformed classical Hermite, Laguerre, Legendre and Jacobi polynomials 
(for the latter, there exists the q-deformati on of the asymmetric form (2.61) only, 
see e.g. Exton [ 1985 ], Gasper and Rahman | 1990[ ), one should modify the spectral 
problem (3.4): 

f k (x,D)<p(x) = eip{qx), (3.8) 

by introducing the r.h.s. function the dependence on the argument qx (cf. (2.5) 
and (3.4)) as it follows from the book by Exton [ 1983fl (see also Gasper and Rah- 
man | 1990| ). Then corresponding q-difference operators having q-deformed classi- 
cal Hermite, Laguerre, Legendre and Jacobi polynomials as eigenfunctions (see the 
equations (5.6.2), (5.5.7.1), (5.7.2.1), (5.8.3) in Exton ]1983[ , respectively) are given 
by the combinations in the generators: 

E 2 = J " J - - {2}J °, (3.9.1) 
E2 = J%j - (f a ' X Jo + (g-°- 1 {a + l})J^, (3.9.2) 
E 2 = -qJ°J° + J- J- + (q- {2}) J °, (3.9.3) 
E 2 = -qo+o- 1 J° J ° + q a J ° J ~ + [q^- 1 - {a}q b - {b}} J ° + {a}J^ , 

(3.9.4) 

respectively. 

Lemma 3.3. If the operator T 2 (for the definition, see (2.15.1)) is such that 
c ++ = and c + = (n — {m})c+o , at some m = 0, 1, 2, . . . 

(3.10) 

then the operator T 2 preserves both V n +i and V m +\, o/nd polynomial solutions in x 
with 8 free parameters occur. 



appears. It reduces the number of independent parameters of the second-order polynomial element 
°f ^ai 2 (R) ? ' I* becomes the standard Casimir operator at q — > 1. 
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As usual in quantum algebras, a rather outstanding situation occurs if the 
deformation parameter q is equal to a primitive root of unity. For instance, the 
following statement holds. 

Lemma 3.4. If a quasi-exactly-solvable operator preserves the space V n +\ and 
the parameter q satisfies to the equation 

q n = 1 , (3.H) 

then the operator Tk preserves an infinite flag of polynomial spaces Vo C V n +i C 

^2(n+l) C • • • C Vk(n+1) C 

It is worth emphasizing that, in the limit as q tends to one, Lemmas 3.1,3.2,3.3 
and Theorem 3.1 coincide with Lemmas 2.1,2.2,2.3 and Theorem 2.1, respectively. 
Thus the case of differential equations in one variable can be treated as a limit- 
ing case of finite-difference ones. Evidently, one can consider the finite-difference 
operators, which are a mixture of generators (3.2) with the same value of n and 
different q's. 
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4. 2 X 2 MATRIX DIFFERENTIAL EQUATIONS ON THE REAL LINE 

4.1. General consideration. This Section is devoted to a description of quasi- 
exactly-solvable , Tfe(x) and exactly-solvable, Efc(i), 2x2 matrix differential op- 
erators acting on space of the two-component spinors with polynomial components 



+l,m+l 



> i • • • ) 

x^x 1 



,x" 



(4.1) 



This space is a natural generalization of the space V n (see (1.5)). The definition of 
quasi- and exactly-solvable operators is also a natural generalization of the Defini- 
tion 2.1 . 

Now let us introduce following two sets of 2 x 2 matrix differential operators: 

T + = x 2 d x — nx + xa + (T~ , 



T° = xd x -2 + l<r + v- 



T~=d x 
n 1 



(4.2) 



named bosonic (even) generators and 

Q 



a 






, Q = 


xa~ 





r+d x 



-a d x 



(4.3) 



named fermionic (odd) generators, where a ± '° are Pauli matrices in standard no- 
tation 



1 






1 



1 

-1 



It is easy to check that these generators form the algebra osp(2, 2): 

[T°,T ± ]=±T ± , [T+,7 1 -] = -2T° , [J,T Q ]=0 ,a = +,-,0 

{Qi,Q 2 } = -T~ , {Q 2 ,Q 1 } = T+ , 
\{{QiM + {Q2M) = +J , \({Qi,Qi}-{Q 2 ,Q2}) = T° , 

{Qi,Qi} = {Q2,Q2} = {Qi,Q 2 } = , 
{Qi,Qi} = {Q 2 ,Q 2 } = (Qi,Q 2 } = o , 



[Qi,T+] = Q 2 , [Q 2 ,T+] = 
[Q 1 ,T+] = , [Q 2 ,T+] = -Q 1 

[Q lt2 ,T°] = ±)-Q h2 



[Q lt T-] = , [Q 2 ,T~] = -Q 1 
[Q 1 ,T"]=Q 2 , [Q 2 ,T"]=0 



[Ql,2,J] 



1,2 ! 



[Ql,2) ^1 — nQl,2 



(4.4) 



This algebra contains the algebra sZ 2 (R) R as sub-algebra. 
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Lemma 4.1. (Turbiner 1994b| ) Consider the space V n +i,n- 

(i) Suppose n > (k — 1). Any quasi-exactly-solvable operator Tk(x), can be 
represented by a kth degree polynomial of the operators (4-2), (4-3). Ifn < (k — I), 
the part of the quasi-exactly-solvable operator Tk(x) containing derivatives in x up 
to order n can be represented by an nth degree polynomial in the generators (4-2), 
(4-3). 

(ii) Conversely, any polynomial in (4-2), (4-3) is a quasi-exactly solvable opera- 
tor. 

(Hi) Among quasi-exactly-solvable operators there exist exactly- solvable operators 
E fc C T fe (x). 

Let us introduce the grading of the bosonic generators (4.2) 

deg(T+) = +1 , deg{J,T°) = , deg(J-) = -1 (4.5) 
and fermionic generators (4.3) 

deg(Q 2 ,Q 1 ) = +\ , degiQ^Qz) = -\ (4.6) 
Hence the grading of monomials of the generators (31), (32) is equal to 
deg{(T+y i +(T ( y a (jf(T-) n -Q 1 mi Q 2 m2 Q 1 ™ 1 Q 2 m2 } = 

(n + — n_) - (ra 1 -ra2-mi+m 2 )/2 (4.7) 

The n's can be arbitrary non-negative integers, while the m's are equal to either 
or 1. The notion of grading allows us to classify the operators Tk(x) in the 
Lie-algebraic sense. 

Lemma 4.2. A quasi-exactly-solvable operator Tfc(cc) C U osp (2,2) has no terms of 
positive grading other than monomials of grading +1/2 containing the generator Q\ 
or Q 2 , iff it is an exactly-solvable operator. 

Take the eigenvalue problem 

T k (x)cp(x) = eip(x) (4.8) 

where 



(p(x) 



.<P2{x). 



(4.9) 



is a two-component spinor. 

Theorem 4.1. (Turbiner ]l994b[ ]j Let n be a non-negative integer. Take the eigen- 
value problem (4-7), where Tk{x) is symmetric. In general, the problem (4-7) has 
(2n + 1) linearly independent eigenf unctions in the form of polynomials V n +x >n , if 
and only ifTk(x) is quasi-exactly-solvable. The problem (4-7) has an infinite se- 
quence of polynomial eigenf unctions, if and only if the operator is exactly-solvable. 

As a consequence of Theorem 4.1, the 2x2 matrix quasi-exactly-solvable differ- 
ential operator T^x), possessing in general (2n + 1) polynomial cigcnfunctions of 
the form V n ^ n -\ can be written in the form 

i—k 

T k (x) = ^T.-u. ,.:r:< , (4.10) 

i=0 
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The coefficient functions &k t i(x) are by 2 x 2 matrices and generically for the fcth or- 
der quasi-exactly-solvable operator their matrix elements are polynomials. Suppose 
that k > 0. Then the matrix elements are given by the following expressions 

/ Ak+i] n[*+»-l] \ 

= ( k D Ui] j (4- 11 ) 

at k > 0, where the superscript in square brackets displays the order of the corre- 
sponding polynomial. 

It is easy to calculate the number of free parameters of a quasi-exactly-solvable 
operator Tfc(x) 

par{T k (x)) =4(fc + l) 2 (4.12) 

For the case of exactly-solvable problems, the matrix elements (4.10) of the 
coefficient functions are modified 

akAx) = [ cip nt ) (4J3) 

where k > 0. An infinite family of orthogonal polynomials as eigenfunctions of Eqs. 
(4.8)-(4.9) , if they exist, will occur, if and only if the coefficient functions have the 
form (4.12). The number of free parameters of an exactly-solvable operator Efe(x) 
and, correspondingly, the maximal number of free parameters of the 2x2 matrix 
orthogonal polynomials in one real variable, is equal to 

par(E k {x)) = 2fc(fc + 3) + 3 . (4.14) 

Thus, the above formulas describe the coefficient functions of matrix differential 
equation (4.8), which can possess polynomials in x as solutions. 

4.2. Quasi-exactly-solvable matrix Schroedinger equations (example). Now 

let us take the quasi-exactly-solvable matrix operator T2 (x) and try to reduce Eq. 
(4.8) to the Schroedinger equation 

h^+V(y)]v]/(y) = E*(y) (4.15) 

where V(y) is a two-by-two hermitian matrix, by making a change of variable x 1— * y 
and "gauge" transformation 

*(y) - U(y)<p(x(y)) (4.16) 

where U is an arbitrary 2x2 matrix depending on the variable y. In order to get 
some "reasonable" Schroedinger equation one should fulfill two requirements: (i) 
the potential V(y) must be hermitian and (ii) the eigenfunctions ^(y) must belong 
to a certain Hilbert space. 

Unlike the case of quasi-exactly-solvablc differential operators in one real vari- 



able (see Gonzalez-Lopez, Kamran and Qlver[ 1993 ), this problem has no complete 
solution so far. Therefore it seems instructive to display a particular non-trivial 
example of the quasi-exactly-solvable 2x2 -matrix Schroedinger operator (Shifman 



and Turbiner |]198J). 

Take the quasi-exactly-solvable operator 

T 2 = -2T°T~ + 2T~J - i/?T°Qi + 
aT° - (2n + 1)T- - y (3n + l)Qi + ^a(3Q 2 - i0Q x , (4.17) 
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where a and (3 are parameters. Upon introducing a new variable y = x 2 and after 
straightforward calculations one finds the following expression for the matrix U in 
Eq. (4.16) 

U = ex P (-^ + ^a 1 ) (4.18) 
and for the potential V in Eq. (4.15) 

V( V ) = i(« 2 - /? V + oa[-(n + \)P +fy*-j tan if] cos ^ + 

^h("+ i )/3+^y 2 - I cot^-]sin^- (4.19) 

It is easy to see that the potential V is hermitian; (2n + 1) eigenfunctions have the 
form of polynomials multiplied by the exponential factor U and they are obviously 
normalizable. 

Recent development of matrix quasi-exactly-solvable differential equations and 



further explicit examples can be found in Brihaye-Kosinski [1994 
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5. Ordinary differential equations with the parity operator 
Let us introduce the parity operator K in the following way 

Kf{x) = f(-x) (5.1) 
The operator K has such properties 

{K,x} + = 0, {K,d x }+ = 0, K 2 = 1. (5.2) 
Now take the following generators 

J~{n) = d x + -(l-K) , (5.3) 

x 

Jl(n) = xJ~ (n) = xd x + v(l - K) , (5.4) 

J°(n) = n — xd x , (5-5) 

J + {n) ~ xJ°(n) — nx — x 2 d x . (5-6) 

These operators toget her wit h the operator K form an algebra, which was named 
gl(2, R)k in Turbiner 1994a] , with commutation relations 

[J°(n),J 1 °(n)]=0, (5.7) 

[J ± (n),J°(n)} = ±Jt(n) , (5.8) 

[J^n), J^n)] = T (l T 2^K) J±(n) , (5.9) 

[J + (n),J-(n)] = Jl{n)-{l + 2vK)J°{n) (5.10) 

and also 

[K,4{n)] = [K,J°(n)] = , {K, J ± (n)} + = 0. (5.11) 

If v = the generators (5.3)-(5.6) and commutation relations (5.7)-(5.10) become 
those of the algebra gl(2, R) []. 

It is easy to check, that the algebra (5.7)-(5.11) possesses two Casimir operators 
- the operators commuting with all generators of the algebra gl(2, ~R)k- the linear 
one 

Ci = J° + J%+vK (5.12) 

and the quadratic one 

C 2 = \{J\ J-} + + \(J° J?) 2 - |(J° - J?)K + \K . (5.13) 

Substituting the concrete representation (5.3)-(5.6) into (5.12)-(5.13) , one can find 
that 

C\ = n + v , 

C 2 = ^[{n + v){n + v + 2) - v 2 ] (5.14) 

The representation (5.3)-(5.6) has an outstanding property: for generic v, when 
n is a non-negative integer number, there appears a finite-dimensional irreducible 



6 Hereafter we will omit the argument n in the generators except for the cases where the 
representation (5.3)— (5.6) is used concretely 
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representation V n +i of dimension (n + 1). So the representation space remains the 
same for any value of parameter v. 

Like it has been done before one can introduce quasi- exactly- and exactly-solvable 
mixed operators, containing the differential operator and the parity operator K : 
Tk(x,d x , K) and Ek(x,d x ,K), respectively. It is evident, that those operators are 
linear in K. 

Lemma 5.1. (i) Suppose n > (k — 1). Any quasi-exactly-solvable operator T k , 
can be represented by a k-th degree polynomial of the operators (5.3) - (5.6) If 
n < (k— 1), the part of the quasi-exactly-solvable operator Ty. containing derivatives 
up to order n can be represented by an nth degree polynomial in the generators ( 5.3)- 
(5.6). 

(ii) Conversely, any polynomial in (5.3)-(5.6) is quasi-exactly solvable. 
(Hi) Among quasi- exactly- solvable operators there exist exactly- solvable operators 
E k CT k . 

Similarly to the case of sZ2(R) one can introduce the grading of generators (5.3)- 
(5.6) 

deg{J + {n)) = +1 , deg(J°(n), J? (n), K) = , deg(J-(n)) = -1, 

(5.15) 

and 

deg[{J + {n)) k +{J°{n)) ko {4{n)) k ^{K) k {J-{n)) k -] = k+ - k- 

(5.16) 

(cf. (2.3)— (2.4)). The grading allows us to classify the operators T k in the algebraic 
sense. 

Lemma 5.2. A quasi-exactly-solvable operator T k C t/ 9 ; 2 (R) K has no terms of pos- 
itive grading, if and only if it is an exactly- solvable operator. 

Theorem 5.1. Let n be a non-negative integer. Take the eigenvalue problem for a 
linear differential operator of the kth order in one variable 

T k (x,d x ,K)ip = 6(p , (5-17) 

where T k is symmetric. The problem (5.17) has (n + 1) linearly independent eigen- 
functions in the form of a polynomial in variable x of order not higher than n, if 
and only ifT k is quasi- exactly- solvable. The problem (5.17) has an infinite sequence 
of polynomial eigen] 'unctions, if and only if the operator is exactly-solvable. If T k 
has no terms of odd grading, (k + — kJ) is odd number (see (5.16)), T k commutes 
with K and eigen) 'unctions in (5.17) have definite parity with respect to x — > —x. 

Following the Lemma 5.1 a general second-order quasi-exactly-solvable differ- 
ential operator is defined by a quadratic polynomial in generators of gl2^R)i<- 
Provided that the conditions (5.12)-(5.14) are taken into account H we arrive at 



7 It leads to a disappearance of the terms containing, for instance, Jj (n) and J°(n)J°(n) 
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T 2 = c ++ J + (n)J+(n) + c +0 J + (n)J°(n) + c+-J + (n)J-(n)+ 

Co-J°(n)J-(ri) + c J-(n)J-(ri) +c + J + {n) +c J°(n) + c_J~(n) + c, 

(5.18.1) 

(cf. (2.15.1), where c Q/3 = C° ap + C^,c a = C° + Cf , c = C° + C K and all C's are 
real numbers. The number of free parameters is par(T 2 ) = 18. Non-existence in T 2 
of the terms of odd grading leads to the conditions 

c+o = c - = 

and 

c + = c_ = 

and, finally, a general operator having eigenfunctions of a definite parity is 

T 2 (e ' o) = c ++ J+{n)J+{n) +c + _J+(n)J-(n) + c__ J-(n)J-(n) + c J°(n) + c, 

(5.18.2) 

and the number of free parameters is pariT^' ^) = 10. 
The condition of Lemma 5.2 requires that 

c++ = c+o = c+ = , 

and then the operator T 2 becomes exactly-solvable 

E 2 = c + _ J + (n)J~(n) +c -J°(n)J _ (n) + c J~(n)J~(n) +c J°(n) + c_J~(n) + c, 

(5.18.3) 

(cf. (2.15.2) and the number of free parameters is reduced to par(E 2 ) = 12. 
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